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1. Introduction 

As a paradigm for the effects of strong, local coulomb interactions, the Anderson 
impurity model (AIM) [1] remains highly topical some forty years after its inception 
(for a comprehensive review, see [2]). Its essential physics in the strong coupling regime 
of large on-site interaction strength (f/) is that of the Kondo effect, characterized by a 
low-energy scale cuk; and manifest famously in the many-body Kondo or Abrikosov-Suhl 
resonance appearing in the single-particle spectrum D{uj). Although ujk itself naturally 
depends upon the interaction strength, the fact that it is the sole low-energy scale means 
that the Kondo resonance exhibits universal (t/-independent) scaling in terms of uj /ojy^ 
alone. 

An obvious question is: what is the form of the universal scaling spectrum? It is 
this we consider in the present paper, in possibly the simplest context of the particle-hole 
symmetric AIM. Perhaps the first point to make is that we do not believe the answer 
to this question is known, and that this reflects in part the well known difficulties in 
constructing approximate theories for dynamical properties of the AIM. The problem is 
of course well understood at low frequencies, where D{ui) — D{Q) oc —{oj /ujy^Y satisfies 
the dictates of Fermi liquid theory, as arises directly from a simple low-frequency 
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expansion of the impurity single-particle Green function [2]. Such behaviour is however 
confined to the lowest of frequencies \uj\/uJ]^ <^ 1. Naive extrapolation of it leads to a 
rather trivial Lorentzian scaling spectrum, as indeed arises in a variety of theoretical 
approaches; for example [2] microscopic Fermi liquid theory to leading order, the slave 
boson mean- field approximation, or D{lj) approximated by the spinon spectrum that 
may be obtained via the Bethe ansatz. 

Numerical approaches by contrast, such as the numerical renormalization group 
(NRG) [3, 4] or quantum Monte Garlo (QMG) [5, 6] calculations, reveal a very different 
behaviour in the form of long, slowly varying tails that entirely dominate the scaling 
spectrum for |a;|/ci;K ^ 1, and are known to be important experimentally [7]. These 
are currently believed [3-6] to be so-called Doniach-Sunjic (DS) [8] tails, of asymptotic 
form D{lj) oc (|c<j|/c<jk)~2 , reflecting physically an incipient orthogonality catastrophe. 
Such behaviour has however been inferred by direct comparison of the numerical results 
to an empirical DS form [3-6]; and while there is no doubt that they are thereby quite 
well described, we do not know of a theoretical approach that (a) exphcitly yields such 
behaviour (if correct), and (b) simultaneously recovers the requisite Fermi liquid form 
as \u!\/u!k — 0. 

We consider these issues within the framework of the local moment approach (LMA) 
that has recently been developed [9-11] to handle in particular dynamical properties of 
AIMs [9, 10]. The aims of the paper are straightforward, and threefold, (i) To obtain 
in closed form the strong coupling LMA scaling spectrum D{uj) for the symmetric 
AIM. This has hitherto been determined numerically in reference [9], and shown to 
give very good agreement with NRG results [4]. But it naturally precluded an explicit 
determination of the form of the dominant spectral tails, which (ii) is our second aim. 
We find that these are not of DS form, but rather exhibit a much more slowly varying 
logarithmic decay. This behaviour is shown to be in large part independent of the details 
of the LMA, and we give further qualitative arguments in support of it. (iii) In view of 
this we reassess comparison between LMA and NRG results [4] for the scaling spectrum, 
concluding in particular that the slow logarithmic tails are very well supported by NRG 
data. We consider in addition the scaling behaviour of the conventional interaction 
self-energy S(c<j), which the LMA correspondingly predicts to diverge logarithmically 
for IcuI/cuk ^ 1- Since recent NRG advances [12] now permit an accurate numerical 
determination of the self-energy itself, this too may be compared to LMA predictions; 
very good agreement is again found. 

§2 of the paper gives a brief introduction to the LMA, and the background required 
for the remainder of the work; full details may be found in references [9, 10]. General 
consideration of the scaling spectrum within the LMA framework is given in §3; and 
in §4 the asymptotic behaviour of the spectral tails is deduced explicitly and compared 
directly to NRG results [4]. In §5 both the LMA scaling spectrum and resultant single 
self-energy are considered on all frequency scales, and likewise compared to results from 
NRG calculations. The paper concludes with a brief summary. 
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2. Background 

The Hamiltonian for the AIM [1] is given in conventional notation by 

-H" = ^ ekflka + XI (fj + ^fli-c^hi^ + X Vik[c\^Cka + h.c.) . (2.1) 



k.a cs k.a 



The first term refers to the host band of non-interacting electrons (with dispersion efe), 
and the second to the impurity with on site interaction U and site-energy e^; for the 
particle- hole (p-h) symmetric AIM considered here, — —U/2 and the impurity charge 

= = 1 for all ^• 

We focus on single-particle dynamics (at T = 0), embodied in the impurity 
Green function G{uj) (<-> G{t) = — i(T{cjo-(t)c|^})) and hence single-particle spectrum 
D{u;) — — 7r~^sgn(a;) Im G{uj). G{uj) is conventionally expressed as 

G(cu) = [cu+- A(a;) -S(a;)]-^ (2.2) 

where uj'^ — u; + iO"'"sgn(u;). Here A(c^') = Ar(cj) — i sgn(cj)Ai(cj) (= —A{—uj)) is 
the host-impurity hybridization, with Ai{u) = 'n'J^k^ik^i^ ~ ^fc)- "^^^ hybridization 
strength Aq = Ai(c<j = 0) oc Phost(^ = 0) is thus defined (with u = the Fermi level), 
and is non-zero since the host is metallic by presumption. S(c<j) = S^(c(j) — i sgn(c(j)S^(c<j) 
is the conventional single self-energy (excluding the trivial Hartree term, which cancels 
ei = -U/2). 

Equation (2.2) merely defines the self-energy S(a;), via the Dyson equation thereby 
implicit. As such it naturally invites a calculation of S(a;) based at heart upon 
perturbation theory in U about the non-interacting limit. But the limitations of such 
approaches - in practice, and certainly in the strong coupling regime - are well known 
(see e.g. [2]); and a determination of G{uj) via the conventional single self-energy is by 
no means mandatory. The LMA [9-11] thus eschews such an approach completely, and 
instead employs a two-self-energy description with G(a;) expressed formally as 

G(u;)^'-[G^(u;) + G^(u;)] (2.3a) 



where 



G,{iu) = [^+ - A(u;) - S,(a;)]-^ (2.36) 



(and a / I or +/—). The interaction self-energies S(t(<^), which by p-h symmetry 
satisfy 

t^(u;)^-t^{-uj), (2.4) 
are separated as 

E,(a;) = -|c/|//|+E,(a;) (2.5) 

into a purely static Fock contribution (with local moment that alone would survive 
at the simple mean-field (MF) level of unrestricted Hartree-Fock; plus an a;-dependent 
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contribution So.(c(;) containing the spin (and charge) dynamics that, at low energies in 
particular, dominate the physics of the problem. 

Use of a two-self-energy description provides a tangible means of developing a 
relatively simple many-body approach to the AIM that starts from, but successfully 
transcends the deficiencies of, the crude MF (or 'frozen spin') approximation. It is 
moreover a necessity and not a luxury for problems that do not ubiquitously exhibit 
Fermi liquid behaviour, but contain an underlying quantum phase transition to e.g. 
a degenerate local moment ground state; the soft-gap AIM [13] provides a specific 
example [10, 4]. The conventional single self-energy S(a;) can of course be obtained as 
a byproduct of the two-self-energy description: direct comparison of equations (2.2,3) 
gives 

H^^iuj) - Et(-a;) + 2g{uj)t^{uj)E^{-u;)} 



E(a;) = 



l-hgiu;)[t^iu)-t^i-u;)] 



(2.6) 



where g{oo) = — A(ci;)]~"^ = g" 
i sgn(a;)7rdo(<^) is the trivial U — propagator. 
From this the gross deficiencies of simple 
MF itself are also seen directly. Here the 
dynamical Scr(c<j)'s are neglected, 



MF, 



S„ = -a/\ 



/-a 



-^U\fi\ (equation (2.5)) and hence Emf(^) 



At the Fermi level in particular 

2 



E[^p(a; = 0) = 7rdo{ui = 0)[^U\iJ,\\" is non-zero, 
thus violating Fermi liquid behaviour, and 
indicative of the broken symmetry inherent at 
MF level (which for |//| 7^ corresponds to a 
doubly degenerate local moment state) . 

The LMA has two essential elements [9, 
10]. First, it includes in the dynamical Eg. (a;) 
a non-perturbative class of diagrams (figure 1) 
that embody dynamical coupling of single- 
particle excitations to low-energy transverse spin fluctuations, and hence capture the 
spin-flip scattering essential to describe the strong coupling/Kondo regime for If — 
U/ttAq » 1. Other classes of diagrams may also be included [9], but retention of the 
dynamical spin-flip scattering processes is essential. These are expressed in terms of the 
MF propagators (sohd lines in figure 1), viz 



Figure 1. Principal contri- 
bution to the LMA S^(a;), 
see text. Wavy lines denote 
U. 



g,{u) = [u;+-A{u;) + ^-U\fi\]- 



[2.7) 



with corresponding spectral densities D^{lj) — — tt ^sgn(a;)Im^o-(<^); and E|(a;) 
E^(a;) — i sgn(a;)Ej(a;) is given explicitly by 



TT 



imu+-{uji)[e{uji)gi{uji+uj) + e{-uji)g+{uji+uj)] (2.8) 
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where 



oo 



Dl{uj,)e{±uj,) 

uj — uji± iO+ 



(2.9) 



are one-sided Hilbert transforms {0{x) is the unit step function). Here. 11+^ (cl^) is the 
transverse spin polarization propagator (shown hatched in figure 1). It is given at the 
simplest level by an RPA-like particle-hole ladder sum in the transverse spin channel; 
viz 



where °n+~(a;) is the bare particle-hole bubble, itself expressed in terms of the MF 
propagators. 

The second, key idea behind the LMA is symmetry restoration: restoration of the 
broken symmetry endemic at MF level, via the spin-flip dynamics embodied in 
This is reflected mathematically in S?(0) = ^^^(0), i.e. (using p-h symmetry) by 



Imposition of equation (2.11) as a self-consistency equation is achieved in practice for 
given U by varying the local moment \^\ from its MF value. It preserves the U- 
independent pinning of the Fermi level spectrum (7rAo-D(c^ = 0) = 1 V f/), and in 
turn leads correctly to Fermi liquid behaviour at low energies [9]. Most importantly 
it introduces naturally a low-energy spin-flip scale cUm - manifest in particular in a 
strong resonance in Imn'^~(a;) centred (by definition of a;ni) on a; = uj-^i - that sets 
the timescale for symmetry restoration. This is the Kondo scale. Its form in strong 
coupling, viz uj-^i oc exp(— ttCZ/SAq) (discussed further below) is asymptotically exact. 

The LMA is readily implemented, as considered in [9] with 11"^ (a;) given by the 
p-h ladder sum equation (2.10). In weak coupling, f/ ^ 0, it is perturbatively exact 
to/including second order in U about the non-interacting limit. More importantly, in 
the strong coupling regime the exponentially narrow Kondo resonance is captured, and 
exhibits universal (^/-independent) scaling in terms of a;/a;m; or, equivalently, in terms 
oi u/uk where the Kondo energy lok (oc tUm) is defined as the HWHM of D(a;). 

The universal scahng regime is reached in practice for [/ > 4 [9], and it is this we 
focus on here. Our aim is to obtain the LMA scaling spectrum analytically, and to do 
so in the first instance with only rather minimal assumptions about the form of the 
transverse spin polarization propagator 11+" (a;). 

3. Scaling spectrum: general considerations. 

To obtain the scaling form of the Kondo/ Abrikosov-Suhl resonance in strong coupling, 
one considers finite uj — uj/uJui in the limit uj^ oc exp(— ttCZ/SAq) — > 0; the Hubbard 
satellites centred on = ^ are naturally not part of the scaling spectrum, and are 



H+-(a;) = °H+-(u;)/[l - C/°H+-(a;)] 



(2.10) 



Ef (0) = -\ 



?7|/i| + Sf (0) = 0. 



(2.11) 
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thereby projected out. Hence, referring to equation (2.3), the 'bare' u = u^^^lo = 
may be neglected, and hkewise A(a;) = A{u!jn^^) reduces to A(0) = — i sgn(a;)Ao. The 
spectrum then follows from equation (2.3) as 

nA„D(.) = lT " + ^°'^^('^)) (3.1) 

2 ^ (A„-'E?(a,))^ + (1 + A„-'Ei(a,))= 

(where S^(a;) = E^(a;)). As expected, its scahng behaviour is determined exclusively by 
that of the interaction self-energies; and we note in passing that equation (3.1) is quite 
general {i.e. provided the host is metallic, it applies for any one-electron hybridization 

A(c.)). 

The LMA "E^^uj) is given by equation (2.8), and in strong coupling the transverse 
spin polarization propagator Imn"'"~(a;) has the following functional form [9] 

-Imn+-(a;) = —f{u;)9{u;) (3.2) 

TT UJ^ 

with 

oo oo 

J ^ Imn+-(a;) = 1 = A^dy f{y). (3.3) 



Three points should be noted here. First, Imn+~(a;) naturally scales in terms of co/ujia- 
Second, equation (3.3) reflects physically the saturation of the local moment (|/x| — > 1) 
and total suppression of double occupancy in strong coupling. Third, the function 
/((D) is peaked at a) = 1 (by definition of cUm), and /(u)) ~ cD as cl; — > 0. The above 
strong coupling behaviour arises explicitly [9] with n+^(c<j) given by the p-h ladder sum 
equation (2.10), from which /(a)) is found to have the form 

/(a)) = -. (3.4) 

In the following, however, we proceed without reference to the specific form of /(cD), 
which will be required only in §5. 

Prom equations (2.8) and (3.2) it follows directly that 

|a;| 

mu;)=9{-^)nU'A J dy f{y)Dl{u;^[y + cI;]) (3.5) 



where D^{uj„i[y + uj]) = D^^iO) since we consider finite uj with uij^ 0. From 
equation (2.7) 7rAo-Dj'(0) = [{^U\n\)'^ + Aq]~^ generally, so in strong coupling {U 1, 

l/^l - 1) 

ttU^DKo) = 4Ao. (3.6) 
Equation (3.5) thus reduces to 

A^'mu) = 9{-u)4A J dy f{y). (3.7) 
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As required from equation (3.1), Aq ^E^(a;) thus scales solely in terms of a) = uj/uj^ with 
no [/-dependence {A is a pure number of order 1, determined via the 'normalization' 
equation (3.3)). 

Likewise, using equations (2.8) and (3.2), 

oo 

Ef (a;) = U'A J dy f{y) Reg;{u;^[y + u]) (3.8) 



the cD-dependence of which is controlled by that of Re^j^(a;) as a; — > 0. Since the latter is 
given from equation (2.9) as a one-sided transform, its low-o; behaviour is logarithmically 
divergent and given by [9] 



X 



\UJ\ 



4Ar 



In 



TT 



A 



\UJ\ 



(3.9) 



(using equation (3.6)). Here, A = min[D, ^] is a high-energy cutoff (with D 
the bandwidth of Ai(a;)); its precise value is immaterial in what follows. From 
equations (3.8,9), 





" A " 








TT 







(3.10) 



The [/-dependence of the Kondo scale ujra now follows from the symmetry 
restoration condition (2.11), viz E|^(a; = 0) = ^ strong coupling. Using 
equation (3.10) this is given simply by 

UJru = C Uj' 



where c is a [/ -independent constant of order unity given by 

" oo 

c = exp -A / dy f{y)\n.{y) 



and 



A exp 



7r[/ 
'8A^ 

■kU 



(3.11a) 



(3.116) 



(3.11c) 



The exponent of oj-^., Viz exp(— ^^), is asymptotically exact in strong coupling (see e.g. 
[2]); the prefactor is of course approximate, and reflects the uv-cutoff used in equation 
(3.9). 

Since T^{u = 0) = (equation (2.11)) has been enforced self-consistently, 
T^{oj) = E^(a;) — E^(0) (see equation (2.5)), and is given from equation (3.10) by 



4A 



TT 



f{y) In 



1 + 



(3.12) 
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As for Aq ^E'j-(c(;) (equation (3.7)), Aq^E^(ci;) also scales solely in terms of d; = uj/u^ 
with no explicit [/-dependence, as required from equation (3.1) for universal scaling of 
7rAo-D(a;) to arise. 

Finally, the quasiparticle weight Z = [1 — ((9S^(a')/(9cj)(^=o]~^) defined in terms 
of the single self-energy E(a;), is also readily obtained. From equations (2.4, 11), 
E^(a;) = E^(a;) oc a; as a; — > 0; and E^u;) (= E^(c<;)) oc u;^ (equation (3.7) using 
f{y) ~ y as y — > 0). Hence from equation (2.6), E^(c<;) = E^(c<;) oc a; as a; — > 0, so 
Z =[1 — {d'E^{u;) /du;)u,=o\~^ is given from equation (3.12) (remembering that LOm — > 0) 

by 

oo 

^"^.^M, (3.13) 



AqZ 7t J y 



Since the problem is characterized by a single low-energy scale, AqZ oc cu^ as expected, 
with the proportionality given explicitly by equation (3.13). 

Equations (3.7,12) (together with p-h symmetry equation (2.4)) provide the 
basic scaling forms for Aq "'^EJ^(c(j) and Aq"'^S^(c<j) that enable the scaling spectrum, 
equation (3.1), to be obtained; they will be evaluated explicitly in §5 for a particular 
form of /(a)). 



4. Spectral tails. 

We consider first the behaviour of the spectral 'tails', = \Lj\/iJ^ ^ 1, since 
the predicted form is not essentially dependent on details of the function /(a)) 
that determines the transverse spin polarization propagator (equation (3.2)). From 
equations (3.7,3), Aq^T.]^-{uj) is given asymptotically for |a;| 1 by 

A^^mu) ^4e{-u). (4.1a) 

Likewise AQ''^E'^(a;) (equation (3.12)) reduces for |a;| » 1 to Ao^E!^(a;) = 
/o°°dy /(y) ln(|a;|/y); and hence using equations (3.3, 11) to 



Aq ^Ef (u;) = — In 



4 



TT 



a; 



(4.16) 



The resultant spectrum for |a;| S> 1 then follows from equation (3.1) as 



7rAoL>(a;) = - i + )> (4.2) 

2\[iin(H)]Vl [iln(M)]V25' 

where cu' — cu/cu'^. 

The scaling spectrum is thus predicted to have a slowly varying logarithmic tail, 
and not Doniach-Sunjic (DS) [8] behaviour of form D{u!) ~ {\uj\/ujk)~'^ ■ The latter has 
hitherto been argued to arise from a numerical renormalization group (NRG) study of 
the asymmetric AIM [3] in the one-hole sector, u; < 0; and, for the symmetric AIM, from 
both quantum Monte Carlo/maximum-entropy studies [5, 6] at finite temperatures and 
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moderate interaction strength, as well as recent (T = 0) NRG work [4] (which data will 
be re-examined below). In all these cases, the tail behaviour was inferred from either 
comparison or fit to an empirical DS form. 

For reasons now explained we believe the logarithmic tails of equation (4.2), 
and not the DS behaviour, to be correct. On physical grounds one expects the 
behaviour of dynamical properties at high frequencies, \u!\/u!k 3> 1, to mirror that 
of static thermodynamic or transport properties at high temperatures, T/Tk ^ 1 
(where the Kondo temperature Tk oc cuk)', and for which asymptotic behaviours are 
characteristically of form [ln(T'/rK)]~"' (see e.g. [2]). A specific relevant example of the 
latter is provided by the impurity resistivity p{T), which at high temperatures probes 
in effect the high-energy tails of the single-particle spectrum. For the antiferromagnetic 
Kondo/s-d model without potential scattering (which for spin S* = ^ is the strong 
coupling limit of the symmetric AIM under a Schricffcr- Wolff transformation), parquet 
resummation leads to the well known Hamann [14] result for p(T)/p(0); its leading high 
temperature behaviour is asymptotically exact and given by 

p(T) _ n'SjS + 1) 1 

p(0) 4 [ln(T/TK)]'' ^ ■ ' 

For the S = j case, equation (4.3) with T replaced by \uj\ recovers precisely the 
asymptotic high-frequency behaviour arising from equation (4.2), viz 7tAoD{uj) = 
D(uj) / D(0) ~ ^[ln(|u;|/u;'jj)]"^ (where to'^^ cx c^k)- It can in fact be shown, starting from 
the general formula [2] for p(r)/p(0) and assuming that the asymptotic high frequency 
behaviour of the single-particle spectrum is given from equation (4.2), that equation (4.3) 
arises directly. 

The second reason we beheve equation (4.2) to be correct is simple: direct 
comparison to the NRG-determined scaling spectrum for the symmetric AIM [4], 
specifically 7iAqD{uj) vs uj/uoy^ (with cjk the HWHM in D{uj)). As discussed in §5 below 
we find for the LMA that c^K/^m ~ 0.691, so that equation (4.2) is directly expressible 
in terms of |cj|/c(Jk- The resultant comparison is shown in figure 2 for uj/uk up to 500. 
The asymptotic behaviour predicted by equation (4.2) gives excellent agreement with 
the NRG data for lu/ujk ^ 10 or so; and we note that the full form of equation (4.2) is 
required for such agreement {i.e. it is not exclusively dominated by its ultimate large-o; 
asymptote ~ [ln(|a;|/a;K)]~^). Since the NRG calculations are of course for finite U 
they naturally show Hubbard satellites at sufficiently large cu/cuk, with deviations from 
universal scaling arising for cu > 0(Ao); this is illustrated in the inset to figure 2 (for 
U ^6). 

We believe the preceding analysis gives sound theoretical grounds for the 
logarithmic tails embodied in equation (4.2), which as above are well supported by 
comparison to NRG calculations. But what of the DS tail behaviour that has hitherto 
been suggested? As shown in reference [4], the NRG spectral tails can indeed be 
empirically well fit to the form nAoD{uj) = a + b{\(x!\/ujK)~'^ ■ Two points should however 
be noted, (i) The resultant fit, while good even up to \u;\/u;k ~ 100, is noticeably poorer 
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Figure 2. Scaling spectrum wAoD{uj) vs w/wk- NRG result (dashed line) 
compared to LMA asymptotic form equation (4.2) (solid line). Inset: on an 
expanded scale out to u/wk = 10** . 

than the form equation (4.2). (ii) A fit to a + 6(|a;|/a;K)~^ requires a 7^ to capture 
the NRG-calculated tail; while the form 6(|c<j|/c(Jk)~2 , which is the strict DS asymptote 
[8], is simply not quantitatively sufficient. Our view in short is that DS behaviour does 
not arise, and that the slower logarithmic decay of equation (4.2) represents the natural 
asymptotic behaviour. 

Finally, note that the asymptotic behaviour of E|(a;) (equation (4.1)) may be used 
in conjunction with equation (2.6) to deduce the corresponding asymptotics of the 
conventional single self-energy this will be considered, and compared to NRG 

results, in the following section. 

5. Scaling spectrum and single self-energy. 

We turn now to the LMA scaling spectrum on all energy scales, to which end we must 
consider the details of /(a)) that determines the transverse spin polarization propagator 
equation (3.2). 

With n+~(a;) given by the p-h ladder sum equation (2.10), /(a;) has been shown [9] 
to have the form equation (3.4), where the [/-dependence of a (and A, see equation (3.2)) 
is given explicitly in reference [9]. From this it is readily shown that in SC U ^ 1, 
a 1 and A ~ [2(1 — a)]^/n 0, such that Af{y) = 6{y— 1) i.e. from equation (3.2) 
;ilmn^ (uj) = 6{u — Ura) rcduccs to a delta function centred on the Kondo scale c^m; 
from now on we refer to this as the LMA(RPA). From equations (3.7,12) it follows 
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directly that 



Ao^Ef(a;) = -Jln|l + a;'| (5.1a) 
Aois;(cj) = 4^(-[l + cu']) (5.16) 

where u' — oj/oj'^^ (and from equation (3.11)). The LMA(RPA) scahng 

spectrum in closed form then follows from equation (3.1); specifically for a;' > (since 
D{uj) = D{-uj)) by 

[|ln|a;'-l|]%[l + 4^(a;'-l)]^/^ 

Equation (5.2) reproduces fully the LMA(RPA) scaling spectrum determined 
numerically in reference [9] (figure 8 therein). For a;' ^ 1 (in practice a;' > 5), it 
naturally recovers the asymptotic behaviour equation (4.2) that is independent of the 
specific LMA(RPA). And the Kondo energy ujy^ is readily determined from the HWHM 
of equation (5.2) to be ojy^/uj'^ — 0.691 (as noted in §4). From equation (3.13) with 
^f{y) = ^{y ~ 1)) = nu^/i, and hence the low-frequency behaviour of the 

LMA(RPA) scaling spectrum follows from equation (5.2) as 

7rAoD{u) ^ — -2- (5.3) 

1 + 

This is precisely the spectrum arising from microscopic Fermi liquid theory to leading 
order [2] {i.e. from the quasiparticle form G{u) ~ [^V-^ + i sgn(cij) Aq]^^). A 
pure Lorentzian arises also from a spinon approximation Ds{uj) to the single-particle 
spectrum, obtained from the Bethe ansatz as [15] Ds{uj)/Ds{0) = [1 + ((^/[^Aq.^^])^]"^. 
However equation (5.3) itself is of course confined to \u;\/AoZ <^ 1 - the Lorentzian 
tails it would otherwise predict are far from correct, as we have seen. 

The deficiencies of the LMA(RPA) D{uj) are not great; as shown in reference [4] 
(figure 10, inset), it gives good agreement with the full NRG scaling spectrum. Its 
principal limitation is for \uj'\ ~ 1 where (see equation (5.2)), a small but anomalous 
'dip' arises in D{uj). This is purely an artifact of the specific RPA-likc form for the 
polarization propagator, reflected in the fact that the resultant lmn+^(ci;) in strong 
coupling becomes a 5-function at a; = ujui' in reality, one additionally expects Imn''"~(a;) 
to have a width on the order of cUm- 

To remedy this deficiency we take a different route from that hitherto considered 
briefly in [9], which will also prove useful in subsequent work on spectral dynamics in 
a magnetic fleld. We retain phenomologically the form equation (3.4) for f{Oj), which 
has both its maximum at Co = uj/uJj^ = 1 for any a G (0,1) and a finite FWHM 
= 2 [(2 — a) ^ — 1] 2 provided a ^ 1; and we employ a high-frequency cutoff uJc (to 
render f{uj) normalizable, equation (3.3)). The parameter a is then determined by 
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requiring that limt^^o(S^(<^)/<^^) is obtained exactly (with S^(a;) the conventional single 
self-energy), i.e. [2] 

(5.4) 



~ 2 



T,(uj) itself is related to S 



IS 



recovered if 



by equation (2.6), and the low-cu behaviour of the latter 
Ao^S}(u;) ~ e{-u)2Au^ (from equation (3.7)), with 
Prom this it is straightforward to show that equation (5.4) is correctly 



-uj/AqZ and Aq^S|(ci;) 



UJr, 



AoZ 



(5.5) 



Equations (3.3, 13) and (5.5) then imply 



fdy ^ 



(5.6) 



which thereby determines a for the chosen cutoff uJc- The latter is of course arbitrary, 
but as expected physically results are not sensitive to it provided it is neither too small 
nor too large. In practice we choose oJc = 10 (a = 0.308); and we note in particular 
that the resultant Kondo energy ljk differs insignificantly from its LMA(RPA) value of 
LOK/'^m ~ 0.691 (with uj'^ from equation (3.11c)). 

The consequent LMA spectrum is shown in figure 3, and is seen to be in excellent 
agreement with the NRG scaling spectrum [4] over essentially the entire frequency range. 
At low frequencies in particular, the exact asymptotic behaviour of D{u!) in strong 
coupling is 



nAoDiu) 



1-5 

2 



[A-^E^(a;)]^ 

2 



.-1^1 



AoZ 



(5.7a) 
(5.76) 



as follows from equation (2.2) using Aq ^E^(a;) ~ —u/AqZ together with equation (5.4) 



for An^E^l 



a; 



This is not fully recovered by the microscopic Fermi liquid form 



equation (5.7) that arises from the LMA(RPA). It is obtained correctly from the above 
LMA (since a is chosen such that equation (5.4) is recovered), although as seen from 
figure 3 (inset) the behaviour equation (5.7) is in practice confined to a very narrow 
frequency domain \u!\/AqZ < 0.1 or so. 
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-10 10 20 30 40 



(Id/ Cd^ 

Figure 3. Scaling spectrum 7rAoI?(a') vs w/wk- NRG result (dashed line) 
compared to LMA (full line). Inset: LMA ttAqD^oj) vs ui/AoZ (solid line) 
compared to exact low-frequency asymptote equation (5.7b) (dotted). 

We also add, as now explained, that equation (5.7) is not recovered in an entirely 
consistent fashion from NRG results [4]. As discussed by Bulla, Hewson and Pruschke 
[12] (BHP), there are two essential ways in which the NRG spectrum may be calculated 
in practice, (i) Directly, with D{lj) obtained as a set of 5-functions with associated 
weights (which are then broadened on a logarithmic scale to recover the continuum); 
this is the traditional approach [16, 17]. (ii) Or, as introduced by BHP [12], the single 
self-energy itself is first calculated as a ratio of two correlation ftinctions. and is then 
used in equation (2.2) to obtain D[u!). This is the method of choice, since in contrast 
to the traditional approach it guarantees the spectral sum rule duj D{uj) = 1, and 
significantly reduces deviations from the Friedel sum rule {7rAoD[uj = 0) = 1) [12]. The 
NRG data shown in figures 2,3 has been thus obtained. 

There are two ways in which the quasi-particle weight Z may then be determined. 
Either directly, from {dT,^{uj)/duj)^j=o = —{Z~^ — 1) (~ —IjZ in strong coupling) 
which defines Z\ or by comparison of AQ^E^(a;) at low frequencies to the exact 
asymptote equation (5.4). The two resultant Z's differ quite significantly in 
strong coupling, typically by some 20% or so for the NRG calculations reported in 
reference [4]. (Alternatively but equivalently, if Z is determined in the natural way 
from (9S^(c<j)/c}u;)tj=0! then equation (5.4) for Aq^S^(c(;) is not correctly obtained.) In 
consequence equation (3.7b) is not recovered consistently, although we emphasize (a) 
that this is of course confined to the low frequency regime, and (b) it should be alleviated 
by more accurate NRG calculations. 
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We now consider the scaling behaviour of the resultant single self-energy S(c<j) that 
arises from the LMA; this follows directly from equation (2.6) once E|(cj) is known. 
The behaviour of S|(c<j) for |c(j| ^ 1 is given by equations (4.1), independently of the 
particular LMA considered. From this one finds the asymptotic high frequency form of 
T,\uj) and S^(a;) to be 



Ao-^S^(^) 



l + -ln^(|u;'|) 



16 



— sgn(^a;J — In \uj'\ 



(5.8a) 
(5.86) 



with uj' — ou/cu'^. The logarithmically divergent behaviour of S(a;) is striking, and (via 
the definition S(a;) = u!~^ — A(cj) — G~^{uj)) is a direct consequence of the asymptotic 
behaviour equation (4.2) of 7rAo-D(cj) that has been shown (figure 2) to give excellent 
agreement with NRG results. Note moreover that equations (5.8) are fully compatible 
with the Hilbcrt transform 7rS^(u;) = J^^dui T}{uji)P{1/{lj — oui)); using which, a 
knowledge solely of S^(a;) for \uj'\ ^ 1 (equation (5.8a)) is readily shown to imply 
precisely the asymptotic behaviour equation (5.8b) for E^(ci;). 




Figure 4. Scaling behaviour of single self-energy, Aq ^I]^(a;) vs uj/lok- LMA result 
(solid line), compared for |a;|/a;K > 10 to high frequency asymptote equation (5.8a) 
(dotted line); and to the NRG result for {/ = 5 (dashed line). 



Figure 4 shows Aq S (cu) arising from the LMA (as specified above with a = 
0.308). The low frequency behaviour Aq^S^(cj) ~ J[u;/Ao2']^ is recovered exactly (by 
construction ), and the asymptotic behaviour equation (5.8a) is rapidly approached for 
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Ic^I/cuk ^ 2-3 or so; the latter is shown exphcitly in figure 4 for uj/ujk > 10 where it 
hes within ~10% of the full ^^^(a;). Figure 4 also shows NRG results for Aq^T,\u;) 
for the particular case U = 5, with which the LMA self-energy is seen to be in very 
good agreement. Deviations of the NRG self-energy from the LMA scaling form set in 
steadily for |a;|/a;K ^ 20-30. This reflects simply the flnite-f/ used, which naturally leads 
to deviations from universahty when uj becomes of order Aq {Aq/uk ~ 80 for U — 5); 
but with progressively increasing U the NRG self-energy falls on the scahng curve out 
to progressively larger values oi cu/cuk- 

6. Summary 

We have considered in this paper a local moment approach to the single-particle 
spectrum D{lj) of the symmetric Anderson impurity model, focussing on the universal 
scaling behaviour characteristic of the strong coupling/ Kondo regime, for which the 
LMA has been shown to provide a simple analytical description. Prom previous 
numerical studies [3-6] it is known in particular that D{u!) contains a long tail that is not 
only an integral part of the scaling spectrum, but in fact dominates its behaviour (the 
crossover to Fermi liquid form occurring only on the lowest energy scales (|u;|/u;k <^ !))• 
The LMA predicts this tail to exhibit a very slow logarithmic decay, rather than the 
power law behaviour hitherto believed [3-6] to arise. This prediction has been shown 
to be very well supported by NRG calculations, and further supporting arguments 
for its form were given. More generally the LMA is found to give good agreement, 
over essentially the entire frequency range, with NRG calculations for both the scaling 
spectrum and conventional single self-energy (the latter in particular being 

predicted to exhibit rather striking logarithmically divergent behaviour for \u!\/u!k ^ 1)- 
The LMA itself is an intrinsically simple, non-perturbative many-body approach [9-11] 
to the dynamics of quantum impurity and related models; the key notion behind it being 
that of 'symmetry restoration' within the framework of an underlying two-self-energy 
description. Despite the unconventional nature of the approach, we believe the results 
presented here provide further evidence for its veracity. 
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